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ABSTRACT 

Aims. Supernova remnants (SNRs) are believed to be the main sources of Galactic cosmic rays. Molecular clouds associated with SNRs can 

produce gamma-ray emission through the interaction of accelerated particles with the concentrated gas. The middle aged SNR W28, for its 

associated system of dense molecular clouds, provides an excellent opportunity to test this hypothesis. 

Methods. We present the AGILE/GRID observations of SNR W28, and compare them with observations at other wavelengths (TeV and '^CO 

(y = 1 — > 0) molecular line emission). 

Results. The gamma-ray flux detected by AGILE from the dominant source associated with W28 is (14 ± 5) x 10"* ph cm"^ s"' for E > 400 MeV. 

This source is positionally well correlated with the TeV emission observed by the HESS telescope. The local variations of the GeV to TeV flux 

ratio suggest a difference between the CR spectra of the north-west and south molecular cloud complexes. A model based on a hadronic-induced 

interaction and diffusion with two molecular clouds at different distances from the W28 shell can explain both the morphological and spectral 

features observed by AGILE in the MeV-GeV energy range and by the HESS telescope in the TeV energy range. The combined set of AGILE and 

H.E.S.S. data strongly support a hadronic model for the gamma-ray production in W28. 
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1. Introduction 

' SNR W28, also known as G6.4-0.1, is a middle aged super- 
. nova remnant (with an age of at least 35 thousand years), lo- 
■ cated at a distance between 1.8 and 3.3 kpc in the inner region 
of the Galaxy {l,b = 6.71, -0.05). W28 is a mixed-morphology 
. SNR, with large (about 48 arcmin) apparent dimensions (see 
' IGreen 20091 1. At radio wavelength a shell structure, where the 
I interaction of the SNR ejecta and the ISM creates a shock, is 
. clearly visible (see Figure [1]). The spectral index between 328 
and 1415 MHz is ff = -0.35 ± 0.18 and shows large local varia- 
tions correlated with the flux density (Dubner et al., 2000 ). The 
'^CO (7 = 1 — > 0) molecular line observations taken by the 
' NANTEN telescope reveals a system of molecular clouds associ- 
ated with the SNR dMizuno & Fukui, 2004). Contours in Figure 
\T\ show the intensity of the ^^CO Une emission integrated over 
the velocity range 3 to 27 km s"' corresponding to the dynamical 
distances compatible with W28. The emission appears concen- 
trated in two main complexes near the north-west (cloud N) and 
south (cloud S) part of the SNR. Evidence of interaction between 
the remnant and the system of molecular clouds is given by the 
detection of 1720 MHz OH maser emission ( |Frail et al., 1994| l 
and by the observation of an unusually high value of the ratio 
CO (7 = 3 - 2)/{J = 2-1) (Arikawa lFai.r 19991 1. 
Two TeV sources have been detected by the HESS Cherenkov 
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telescope near W28 ( |Aharonian et al., 2006| l. The sources HESS 
J1801-233 and HESS J1800-240 (A, B and C) are very weU 
correlated with, respectively, the N and S molecular clouds re- 
solved by the NANTEN telescope ( Aharonian et al. , 2008| ). In 
the MeV-GeV energy range a source compatible with W28 
was reported in the Third EGRET catalog, 3EG J 1800-2338 
( |Hartman et al., 19^9| ). The flux measured by EGRET above 
100 MeV was 61.3 + 6.7 ph cm"^ s"', with a photon spec- 
tral index of 2.10. In both the Fermi LAT Bright Source 
List (OFGL J1801.6-2327, see |Abdo et al., 2009|l and the Fi rst 
AGILE Catalog (lAGL J 1803-2258, see |Pittori et al, 2009| l a 
gamma-ray source centered on the position of cloud N is re- 
ported. 

We present in this Letter the results of deep AGILE observa- 
tion of lAGL J1803-2258 which, at energies greater than 400 
MeV, has a shape that turns out to be incompatible with a single 
point source, and is remarkably coiTelated with both the molecu- 
lar cloud system seen by NANTEN and the TeV sources detected 
by HESS. 



2. Data Analysis 

The Gamma-Ray Image Detector (GRID) on board AGILE 
satellite ^Tavani et al., 2 008^ extensively observed the Galactic 
plane during the years 2007-2009 at energies greater than 100 
MeV with an angular resolution better than the previous gamma- 
ray telescopes. Level- 1 AGILE-GRID data were analyzed us- 
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Fig. 1. Map of VLA 90 cm radio emission from SNR W28 in 
celestial J2000 coordinates R.A. and Dec. (Colours indicate the 
intensity in Jy beam '). The black contours show the CO inten- 
sity emission obtained with the NANTEN radio telescope which 
trace molecular clouds, integrated over the velocity interval of 
3-27 km s"' (contours vary between 90 and 170 K deg with a 
step size of 10 K deg). Two concentrations of CO emission are 
clearly visible near the positions a,S = 270.4, -23.4 (cloud N) 
and a, 6 = 270.2, -24.1 (cloud S). 



ing the AGILE Standard Analysis Pipeline. The Anticoincidence 
System (ACS) and a set of hardware triggers perform the first 
reduction of the high rate of background events (charged par- 
ticles and albedo gamma-rays) interacting with the instrument. 
A dedicated software processes the GRID and ACS signals and 
performs a reconstruction and selection of events in order to 
discriminate between background events and gamma-rays, de- 
riving for the latter the energy and direction of the incoming 
photons. A simplified version of this software operates directly 
on board AGILE in order to reduce the telemetry throughput 
( [Giuliani et al., 2006) 1. A more complex version of the photon 
reconstruction and selection software is applied on the ground, 
producing a photon list containing arrival time, energy and di- 
rection of every gamma-ray and the corresponding exposure 
history. Counts, exposure and Galactic background gamma-ray 
maps, the latter based on the Galactic diffuse emission model 
developed for AGILE ( [Giuliani et al., 2004) , were created with 
a bin-size of 0.°05 x 0.°05 for photons with energy greater than 
100 MeV and 400 MeV. We selected only events flagged as 
confirmed gamma-ray events. All events collected during the 
South Atlantic Anomaly or whose reconstructed directions form 
angles with the satellite-Earth vector smaller than 90° are re- 
jected. In order to derive the average flux and spectrum of the 
source we ran the AGILE point source analysis software ALIKE 
( [Chen et al., 2010| ), based on the maximum likelihood technique 
as described in Mattox et al. (1993), over the whole observing 
period. All the flux errors found by ALIKE (and reported in the 
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Fig. 2. Sky map in celestial J2000 coordinates, R.A. and Dec, of 
the significance (cr) of the TeV emission detected by the HESS 
telescope ( [Aharonian et al., 20081 . The blue circle indicates the 
location of the supernova remnant W28. As in Figure[T]the black 
contours show the CO intensity emission. 



following) take into account only the statistical uncertainities. 
We estimate that the systematic errors are on the level of 10% of 
the reported fluxes. 

2.1. Results 

A prominent gamma-ray source associable with W28 (lAGL 
Jl 803-2258) is clearly detected at a significance level of 7 cr 
with an average flux of (40 + 1 1 ) x 1 0"^ ph cm"^ s" ' for E > 1 00 
MeV above the prediction of the AGILE diffuse emission model 
. The positional analysis gives an elliptical error-box centered 
in / =6.0° b =-0.2° with a mean radius of about 0.1° consistent 
with both the EGRET and Fermi sources detected in the same 
region. 

Figure[3]shows the counts map for the source lAGL J 1803-2258 
for energies greater than 400 MeV. At these energies the good 
angular resolution of the GRID detector allows to perform a mor- 
phological analysis of the source, which shows a remarkable cor- 
respondence with the TeV emission observed in the same region. 
The total flux of the source at energy greater than 400 MeV is ( 14 
+ 5) X 10"^ ph cm"^ s"'. Most of the emission is concentrated 
in a region coincident with the molecular cloud N (and the TeV 
source HESS J1801-233). A weaker gamma-ray diffuse emis- 
sion appears to be superposed on the molecular cloud S (and the 
TeV source complex HESS J 1800-240). Performing a likelihood 
analysis with two point sources fixed at the locations of source 
N and source S, we find flux values F^ = (10 ± 3) x 10"** ph 
cm"^ s"' and Fs = (4 ± 2) x 10"** ph cm"^ s"' above 400 MeV. 
The two regions are not clearly resolved in the 100-400 MeV 
energy range. Nevertheless, we performed a likelihood analysis 
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assuming the same two point sources as in the analysis for the 
E > 400 MeV energy range. This yields a flux of (30 + 6) x 
10^^ ph cm"^ s ' for cloud N and a 2cr upper limit of 10 x 10 ** 
ph cm"^ s ' for cloud S. 
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Fig. 3. Gaussian-smoothed AGILE counts map in Galactic co- 
ordinates for the W28 region (counts bin '). Only photons with 
energy greater than 400 MeV have been used. The blu circle in- 
dicates the radio location of the supernova remnant W28. As in 
Figure[T|the black contours give the CO intensity emission. 



3. Discussion 

Figure |4] shows the spectra for cloud N (black) and cloud S (red) 
in the GRID and in the HESS energy bands. While in the TeV 
energy band the complex of the source HESS 1800-240 is a fac- 
tor of two brighter than the source HESS 1801-233, in the GeV 
band the reverse is true. Within a hadronic scenario, and assum- 
ing CRs are accelerated by W28 and then escaping, the differ- 
ence between the TeV to GeV flux ratios determined at the loca- 
tions N and S implies diff'erent proton spectra in the two clouds. 
This is not surprising because, for a turbulent medium, the diffii- 
sion coefficient of charged particles is expected to depend on the 
particle kinetic energy, D = D(E). As a consequence, in the pres- 
ence of recent acceleration (less than 10^ years), the proton spec- 
trum depends strongly on the distance from the accelerator itself 
|Aharoman & Atoyan, 19961 |Gabici et al, 2009| l. Hence the dif- 
lerence in the gamma-ray spectra of cloud N and S can be ex- 
plained by assuming different distances between the clouds and 
the SNR. In the following, we show that this model can explain 
the GeV/TeV spectra of the clouds assuming a set of parameters 
(shown in Table 1 ) well compatible with other observational con- 
straints. We estimated (using the mathematical formalism devel- 
oped in Aharonian & Atoyan, 1996 1 the evolution of the protons 
and nuclei spectra diffusing in the interstellar medium at differ- 



ent distances from the supernova remnant. We assume a distance 
of 2 kpc and an age of the SNR of 45000 years and that the parti- 
cles have been injected continuously during the first 10000 years 
of its life, with a power-law energy spectrum of index 2.2. After 
the injection we assume that they propagate in the interstellar 
medium with diffusion coefficient given by : 
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with 6 - 0.5 and a rather slow diffusion coefficient Do - 10 
cm- s"' as indicated by Aharonian and Atoyan (1996) for a 
dense medium such as that of W28. 
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Fig. 4. Gamma-ray photon spectrum for W28 (AGILE points: 
Icr statistical errors, HESS points are from Aharonian et al. 
2008). The curves represent the gamma-ray spectra estimated, 
accordling to the model presented in the text, for the cloud N 
(black) and the cloud S (red) assumed at a distance of 4 and 9 
pc, respectively, from the particle acceleration site. 



In Figure|4]we compare the results of this model with the AGILE 
and HESS spectral data. The dashed lines correspond to a cloud 
accelerator distance of 4 (in black) and 9 (in red) parsecs. Note 
that a different set of parameters, the cloud-distances and the 
SNR-age, can lead to similar results (in general, any set of these 
parameters with the same ratio distances/ ^JSNR age can repro- 
duce the curves shown in Fig|4|. We choose here the distance of 
the cloud S as the minimum distance compatible with the obser- 
vations. A distance of 4 pc for cloud N (which requires assum- 
ing a projection effect of a few parsecs) and an age of 45,000 
years for W28 are then required in order to reproduce the GeV 
to TeV ratio of the clouds. We set the mass of the cloud S to be 
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larger (by a factor 2.25) than that of cloud N, as can be inferred 
from the TeV fluxes and from the proton density predicted by 
this model at 4 and 9 pc from the SNR shell. The distances used 
in this model are consistent with the position of the clouds as 
observed in the radio maps (assuming an Earth-SNR distance of 
2 kpc, 10 parsecs corresponds to 0.29 deg), taking into account 
posible projection eff'ects. Assuming a mass for the two molec- 
ular clouds, it is possible to evaluate the total energy carried by 
the protons accelerated by W28. We used 2.0 x 10'*Mq for cloud 
N and 4.5 x 10"* M0 for cloud S, in agreement with the estimate 
based on NANTEN data given in [Aharonian et al., 2008] This 
leads to a proton total energy of 3.3 x 10'*'' erg, a value in good 
agreement with the CR energy release from a SNR expected by 
theoretical considerations. 

Finally, we briefly consider here a leptonic scenario for the 
gamma-ray emission from W28. Relativistic electrons can 
produce gamma rays through inverse Compton (IC) and/or 
Bremsstrahlung in a dense target. The cooling times for these 
processes are : 



ttr =4X10^ 



1 cm 



yr 
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4. Conclusions 



The AGILE deep observations of SNR W28 determine the exis- 
tence of several contributions to the gamma-ray emission from 
this SNR. A highly significant gamma-ray source (7 cr for E> 
100 MeV) is associated with W28. This source has an average 
flux of (40 + 1 1) X 10 ** ph cm"^ s"' during the period covered 
by the AGILE observations. AGILE also detected, for photon en- 
ergies above 400 MeV, a SE extension corresponding to a mas- 
sive molecular cloud. We proposed a model based on a hadronic- 
induced interaction with two molecular clouds adjacent to the 
SNR in order to fit the observations. This model explains the 
morphological and spectral features detected both by AGILE i 
n the MeV-GeV energy range and by the HESS telescope in the 
TeV energy range. Setting the distances and masses of the two 
main molecular clouds to values compatible with the radio CO 
observations we can estimate that the total energy in protons is 
3.3 X 10"*^ erg. Our data and analysis provides strong support to 
a hadronic origin of the gamma-ray emission from W28. 



t,c(E) = 3 X lO** 
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where n is the gas density and w,- the energy density of the radia- 
tion field (note that eq. 3 refers to IC scattering in the Thompson 
regime). For comparison the cooling time for protons producing 
gamma-rays by proton-proton interactions is ; 



tpp = 1.6x10° 



1 cm 



yr 



(4) 



Inverse Compton is disfavoured becausethe strong detected cor- 
relation (especially at TeV energies) between the gamma ray in- 
tensity and dense gas is not expected for evolved SNR like W28. 
Bremsstrahlung can, in principle, produce the observed gamma- 
ray emission. Assuming that TeV gamma-rays are produced by 
~ 10 TeV electrons, Bremsstrahlung dominates over the IC in 
very dense regions (« > 1000 cm"^). However such energetic 
electrons are expected to be produced only at the early epochs 
of the SNR ( |Gabici et al., 2 009), and they can only marginally 
survive for the age of this source due to strong synchrotron and 
IC cooling. Furthermore, in order to explain the observed fluxes, 
a leptonic model would require a very large e/p ratio in acceler- 
ated particles (larger than 0.25, comparing eq. 2 and 4) while a 
ratio near 0.01 is expected from the CR measurements. 

Table 1. Parameters of the model used to fit the AGILE and 
HESS spectral data, see figlD 



SNR age 


45000 yr 


SNR distance 


2 kpc 


CR total energy 


3.3 X lO-*' erg 


CR Injection index 


2.2 


Diffusion coeff. @ 10 GeV (Do) 


10^* cm^s-' 


Diffusion coeff. spectral index (d) 


0.5 


Cloud N (mass, distance from the W28 shell) 


2.0 X lO*M0,4pc 


Cloud S (mass, distance from the W28 shell) 


4.5 X 10*Mo,9pc 
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